Universal Heat Transport in Sr 2 Ru04 



(N 
O 

o 

< 



c 

o 
o 

I 

5-H 

Or 

3 : 

. 



I 

C 

o 
o 



> 
o 



o 

(N 

o 

-I— > 



i 

C 

O 

o 



M. Suzuki 1 , M. A. Tanatar 1 ' 2 !, N. Kikugawa 1 ' 2 , Z. Q. Mao 1 ' 2 , Y. Maeno^ 2 ' 3 0, and T. Ishiguro 1 ' 2 
1 Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
2 CREST, Japan Science and Technology Corporation, Kawaguchi, Saitama 332-0012, Japan 
3 International Innovation Center, Kyoto University, Kyoto 606-8501, Japan 

(February 1, 2008) 

We present the temperature dependence of the thermal conductivity k(T) of the unconventional 
superconductor S^RuCU down to low temperatures (~100 mK). In T — ► K limit we found a finite 
residual term in k/T, providing clear evidence for the superconducting state with an unconventional 
pairing. The residual term remains unchanged for samples with different T c , demonstrating the 
universal character of heat transport in this spin-triplet superconductor. The low-temperature 
behavior of k suggests the strong impurity scattering with a phase shift close to tt/2. A criterion 
for the observation of universality is experimentally deduced. 



PACS numbers: 74.20.Rp, 74.25.Fy, 74.70.Pq 

Through the study of the cuprate superconductor in 
the last decade, a notable progress has been made in 
an understanding of the thermodynamic and transport 
properties of quasi-particles (QP) in unconventional su- 
perconductors. One of the important findings is the 
novel phenomenon called the universal transport, pointed 
out first by Lee In an unconventional superconduc- 
tor with nodes in the superconducting (SC) gap, non- 
magnetic impurities suppress the transition temperature 
(T c ) and induce the finite density of QP states at the 
Fermi level with the energy width 7i mp . The presence 
of impurities increases both the QP density and the impu- 
rity scattering rate, which completely cancel each other 
HI}. In this case the QP transport, such as thermal con- 
ductivity (k), restores the temperature (T) dependence 
typical for the normal state (e.g. k oc T) § and be- 
comes independent of impurity concentration. The mag- 
nitude of the residual term in k/T depends only on the 
QP spectrum near the gap nodes, providing a useful tool 
for studying the SC order parameter [Q . 

Experimentally, the universal residual term in k/T at 
T — ► K has been reported in optimally-doped high- 
T c cuprates YBa 2 Cu 3 6 . 9 @ and Bi 2 Sr 2 CaCu 2 8 
Their residual values are consistent with the d-wave SC 
gap structure ||. For non-d-wave gap symmetries, the 
universal conductivity has never been observed. In con- 
trast, the extrapolated k/T at T — > K of the spin-triplet 
superconductor UPt3 [nOUO] rapidly increases with the 
density of defects, showing no universal behavior, and 
the residual term seems to vanish in zero-disorder limit 
Ji"i"| |. The origin of the lack of universality is still unclear, 
similar to the case of under-doped cuprates j^J , and it is 
strongly desirable to extend studies to other unconven- 
tional superconductors. 

In this respect, of special interest is the study of 
Sr 2 Ru04, a unique quasi-two-dimensional (Q2D) spin 

This material has a lay- 
and its electronic state is 



triplet superconductor |1 
ered perovskite structure 
well described in the Landau- Fermi liquid model Jl6| . 
The SC state of Sr 2 Ru0 4 (T c = 1.5 K) is chiral with a 




broken time reversal symmetry [ jL7] and its T c is rapidly 
suppressed by non-magnetic impurities [ p~8|JTs| ] . Views on 
the SC gap symmetry of Sr 2 Ru04 changed in time. Early 
theoretical prediction of the fully gapped p-wave chiral 
state po|| did not find support in recent experiments, sug- 
gesting the SC gap with line nodes pi] p4[ , likely running 
parallel to the conducting plane 1 25 ^27]] . Several new the- 
oretical models have been proposed p8|-|3l]] to explain the 
experimentally observed " triplet- cbiral-nodal" SC order 
parameter. 

The thermal conductivity measurement of Sr 2 Ru04 
down to 20 mK was reported by Suderow et al. J32| . This 
study was restricted to strongly disordered samples with 
T c < 1 K. Measurements with high quality samples, with 
T c > 1.4 K, have been performed only down to 0.3 K 
pd] pfj[ , which is insufficient for the estimation of the 
residual term. 

In this Letter we report the systematic study of the in- 
plane thermal conductivity of Sr 2 RuC>4 down to 100 mK 
on samples with widely varying quality, with T c from 0.7 
K to 1.5 K. We found a universal behavior of the thermal 
conductivity in this material with non-ci-wave order pa- 
rameter. The finite and universal k/T at T —> K in the 
zero impurity limit provides a key support for the SC gap 
with line nodes. Negligibly small phonon contribution to 
the total k allowed us to characterize the T dependence 
of the electronic k, which was not possible up to now in 
the high-T c cuprates ||. We found a reasonable agree- 
ment with the theoretical predictions by Graf et al. 0J in 
the unitarity scattering limit. 

Single crystals of Sr 2 Ru04 were grown in an infrared 
image furnace by the floating zone method (3^] . We used 
samples of 5 different batches with T c = 1.44 K, 1.32 K, 
1.27 K, 1.09 K and 0.71 K, denoted as # 1 - 5, respec- 
tively. The T c and the transition width ST C were defined 
at the peak temperature and the full width at half max- 
imum of the imaginary part of the AC susceptibility, re- 
spectively (Table I). The difference in T c comes from the 
variation of the density of non-magnetic impurities, such 
as Al or Si fill, and/or crystalline defects ]19j . The sys- 
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tematic studies on the electrical resistivity p and resistive 
T c revealed that the suppression of T c is well described 
by the Abrikosov-Gorkov type equation Q, 

M^) = n\ + ^Tjr) ~ Hi), (1) 
l c 2 Airhole 2 

where ^(x) is the digamma function, T c o = 1.5 K is max- 
imum T c for the disorder free material and T is the non- 
magnetic impurity scattering rate in the normal state. 
The values of T for the samples under study, deduced 
from Eq. ([!]), are listed in Table I. 

The thermal conductivity was measured along the in- 
plane [100] direction between 100 mK and 2 K by the one- 
heater and two-thermometer method described elsewhere 
. Typical dimensions of the samples were 2.0 mm x 0.5 
mmx0.08 mm, with the longest direction corresponding 
to the [100] crystal axis and the shortest to the [001] axis. 

Fig. [j] shows the T dependence of k/T of samples # 
1 to 5, with the T c indicated by arrows. All the sam- 
ples show almost flat T dependence of k/T in the normal 
state, as expected for the impurity scattering regime, ex- 
cept for sample # 1 showing slight decrease with T in- 
dicative of small electron-electron scattering. The normal 
state conductivity increases with T c , i.e. with decreasing 
r. Below T c , k/T decreases monotonically and slowly. 
In sample # 1, k/T at 0.3 K (0.2 T c ) is - 5 W/K 2 m (0.2 
k(T c )/T c ), which should be practically zero for the fully 
gapped superconductor Ipq] . Although the normal state 
k/T decreases 4 times with increasing T, in the SC state 
all k/T curves converge to a finite value at around 2 ± 
0.5 W/K 2 m at very low T. 

The p and k/T measurements made on the same con- 
tacts in zero H above T c revealed that the Wiedemann- 
Franz ratio pk/T at T c is satisfied within accuracy of our 
determination (±10% Q) for all sample qualities under 
study, showing a dominance of the electronic transport 
in the normal state at temperatures as high as T c . In the 
SC state, the phonon conductivity re g can be estimated in 
the boundary scattering limit, providing an upper bound 
for k s , by using a simple kinetic formula k^ = ^C p hvl, 
where C p h is the phonon specific heat per unit volume 
proportional to T 3 , v is the averaged sound velocity and 
I is the phonon mean free path, which is of the order of 
the sample dimension. Since both v and I are tempera- 
ture independent, it can be expressed as k^/T = AfT 2 
with a constant coefficient Af . Taking the phonon spe- 
cific heat coefficient of 0.197 mJ/K 4 mole [|l] and w=1400 
m/s [0@1> we obtain Af = 0.2 W/K 4 m, which is a small 
fraction of the total k in all samples. 

On T increase, k/T vs. T deviates from the univer- 
sal value, shows an upward curvature up to 0.3 K and 
then turns into a linear increase at higher T. The k/T 
vs. T 2 dependence (Fig. |l| inset) shows that super-linear 
portion of the k/T vs. T curve is well reproduced by a 
polynomial expression k/T = + X^T 2 (we show only 
the data for samples # 1 and # 3 for clarity) . We should 
notice here that, despite having the same T 2 dependence 



as expected for k 6 /T, the observed k/T is completely 
electronic in origin, which is of special interest in this 
study. Actually measured A| decreases with decreasing 
T c from 44 W/K 4 m for sample # 1 to 8 W/K 4 m for 
sample # 5, but still is 40 times larger than Af of 0.2 
W/K 4 m. The finite residual coefficient Aq indicates ex- 
istence of QP at T=0 and, consequently, the non-s-wave 
pairing. In discussion of the residual term which follows 
below, we use a conventional notation, Aq = kqq/T. 

In Fig. H we plot (filled squares) k q /T as a function of 
the impurity scattering rate divided by the impurity-free 
T c , hr/ksTcQ. The vertical error bars reflect the scatter 
of observed k/T and the horizontal error comes from the 
transition width ST C . For the reference, we show with 
filled triangles k/T by Suderow et al. (3^] extrapolated at 
T — > K. It can be clearly seen that the residual thermal 
conductivity Koo/T remains unchanged for T%T/k B T c o < 
0.2. Moreover, Koo/T values converge to a finite value 
of 1.7 ±0.15 W/K 2 m in the hr/k B T c0 -> limit, giving 
strong evidence for the nodal gap structure. The finite 
kqq/T should appear even in the superconductors with 
fully gapped non-s-wave state (e.g., isotropic 2D p-wave 
state p0[ ]) when a small number of impurities exist, but 
it should vanish in the T — > limit with breaking of the 
universality fl39| . 

Theoretically, the universal residual term Koo /T for the 
nodal superconductor is explicitly expressed as 

where Nf and v-p are the density of states in the nor- 
mal state and the Fermi velocity, respectively. Ao is the 
maximum amplitude of the gap, p is the slope of the 
gap at the node defined as p = > where (j) is the 

in-plane angle of the wave vector on the circular Fermi 
surface, and a is the coefficient of the order of unity which 
depends on the gap symmetry Q. Thus from the resid- 
ual value koo /T we can obtain the quantity related to the 
gap structure, juAo/ci. 

By assuming that the SC gaps have the same magni- 
tude on the three Fermi surfaces of Si^RuGi, and tak- 
ing experimentally measured Fermi velocities and effec- 
tive masses for respective bands Eq ], we can reproduce 
the observed residual value of Koo/T — 1-7 W/K 2 m with 
pA /a = 0.39 meV. This value is close to typical values 
for several nodal gap structures. For the 2D /-wave gap 
with four vertical line nodes A = Ao cos(20) |30[] , a is 4/7T 
M, p = 2, Ao = 2.14/jbTc in the weak coupling limit and 
we obtain p,Ao/a = 0.43 meV. For the state with hor- 
izontal line node, indicated in recent directional exper- 
iments [p5|-p7|, the universal thermal conductivity was 
not considered theoretically. However, we can estimate 
the residual value of koo /T by assuming that in the limit 
of the small Fermi velocity along the [001] axis vp± — > 
0, the 3D polar state A = Ao cos(8) (9 is polar angle 
measured from the [001] axis) is a good representation of 
a Q2D horizontal nodal state. By using the parameters 
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for the polar state | a=l, /i = 1 and A = 2A6k B T c , 
we obtain /uAo/a = 0.32 meV, which is also close to the 
experimentally observed value. Thus we conclude that 
the observed residual term comes from the contribution 
around the gap nodes. 

Recently Zhitomirsky and Rice proposed the or- 
bital dependent gap structure, which is a modified exten- 
sion of the original theory by Agterberg et al. [^8| . They 
assumed the line nodes only on two of three Fermi sur- 
faces, a and (3. If we consider that only these two sheets 
make contribution to k 00 /T — 1.7 W/K 2 m, we obtain 
/iAo/a = 0.27 meV. They also assumed a very small sub- 
Kelvin amplitude of the gap A on these sheets; thus to 
explain the observed value of Koo/T we need a relatively 
steep slope n at the nodes. 

Since the observed k/T is almost purely electronic, we 
can compare the experimental results with available theo- 
ries for the electronic thermal conductivity. Two theories 
predict a T 2 term in the electronic k/T. The model by 
Zhitomirsky and Walker |^0|and by Graf and Balatsky 
|d| (ZWGB) is valid in the range of 7 imp < k B T due to 
the QP density coming from the nodal SC gap in case of 
unitary scattering. The model by Graf et al. Q considers 
finite T correction to the residual conductivity koo /T in 
the range of k B T < jimp, where the coefficient of the T 2 
term of k/T strongly depends on the impurity scatter- 
ing phase shift 0. As can be seem from Inset in Fig. 
[j], the T 2 term is observed from the lowest T, directly 
above universality limit. This implies that the condition 
7imp < k B T required for observation of the T 2 term in 
ZWGB model is not met. Therefore, we stick in inter- 
pretation to the Graf model Q| . Here the T 2 term is ob- 
served as long as we satisfy the condition for k B T < 7i mp , 
which corresponds to T < 0.3 K in our experiment. This 
relatively large impurity band means that the scatter- 
ing is strong. Indeed, in the unitarity limit (the phase 
shift 5 = 7r/2) the impurity bandwidth increases rapidly 
with the impurity scattering rate as 7i mp ~ y/AohT [Q, 
which corresponds to about 0.5 K for sample # 1, while in 
the Born limit (5 — 0) the impurity bandwidth becomes 
exponentially small, 7i mp ~ Aoexp(— ^) (|]. Experi- 
mentally, the impurity bandwidth can be estimated from 
the parameters of the fitting formula, k/T — Ag + A|T 2 
introduced above. In the unitarity scattering limit, the 
deviation of k/T from the universal value at finite T is 
described as B 



K 

T 



K 00 

T 



(1 



i£<— )2) - 

60 7 imp 



(3) 



The 7im P deduced by comparing Eq. (g) with the coef- 
ficients Ag and A§ for each sample is plotted in Fig. |2| 
with an open circle. The uncertainty due to the possible 
small phonon contribution is at most 10 %. The least 
square fit with the %/T 7 function 7i mp = 0.7 \f%T /k B T c0 
is drawn with a dotted line. The rapid increase and \/T 
dependence of 7i mp support the strong impurity scatter- 
ing and verify the deduction of 7i mp from Eq. (||), while 



the absolute value is about two times smaller than the 
theoretical value of 7i mp ~ V AqHT. A slight deviation 
of S from 7r/2 makes the impurity bandwidth exponen- 
tially small at small T B and thus the non-negligible 
7i mp of sample # 1 (UT /k B T c0 = 0.051) implies the large 
scattering phase shift <5 = 7r/2(±10%). 

From our data we can derive the explicit condition of 
the low T and small 7i m p "universality limit" which is 
theoretically described as k B T < 7; mp « A . In Fig. ^ 
we can see that above TiT /k B T C Q — 0.2 the universality 
starts to break and this crossover point corresponds to 
7im P /fcB ~ 0.3 K ~ 0.lA o /fc B . Thus k B T < j imp < 
0.1 A can be considered as an experimentally deduced 
criterion for the observation of universal conductivity. 

The deviation of kqq/T from the universal value 
with increasing HT/k B T C Q is small, even for the largest 
hT/k B T c0 , giving at most 50 % increase for sample # 5 
(TiT /k B T c o — 0.6). According to the calculation by Sun 
and Maki |3| in the unitarity scattering limit, Koo/T in- 
creases rapidly with T above the universality limit and it 
reaches almost twice of the universal value at TiT /k B T C Q 
= 0.6. Similar slow increase of kqq/T was also pointed 
out in YBa2Cu306.g |§- While the slight deviation of the 
scattering phase shift from 7r/2 might be able to explain 
it Q, the calculated steep increase of Koq/T || is inap- 
propriate for the description of Si^RuO,!. The critical 
resistivity, above which T c disappears, is about 1 /ificm 
fl8| , [l9f and the corresponding critical k/T is obtained by 
the Wiedemann-Franz law as 2.4 W/K 2 m, which is al- 
most equal to our largest kqq/T. Since there is no clear 
reason for k/T at T=0 K to become higher in the SC state 
than in the normal state, the increase of kqq/T with T 
expected in Ref. || should be limited. This point should 
be clarified in the further studies. 

In conclusion, we experimentally observed universal 
thermal conductivity in S^RuGi, indicating that this 
phenomenon is not restricted to c?-wave gap symmetries. 
The universality gives clear evidence for the nodal gap. 
From the analysis of the low temperature thermal con- 
ductivity within the framework of the theory by Graf et 
al. [Q, we extract the impurity bandwidth and the scat- 
tering phase shift close to ir/2. We experimentally derive 
a criterion for the deviation from the universal limit. 
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TABLE I. The properties of 5 single crystals of Sr2Ru04. 
T c and ST C were determined from the AC susceptibility mea- 
surements. The impurity scattering rate normalized by the 
maximum T c , hr/ksTco is deduced from Eq. (|l|). 
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FIG. 1. The temperature dependence of the thermal con- 
ductivity k/T along the [100] direction for 5 samples of differ- 
ent quality. The arrows denote the transition temperatures 
T c . Inset shows the k/T vs. T 2 plot for samples # 1 and # 
3. 
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FIG. 2. The residual thermal conductivity Koo/T and the 
impurity bandwidth 7i mp as a function of scattering rate 
HT /ksTco. Filled squares show Koo/T of this study, filled tri- 
angles are by Suderow et al [32]. The open circles represent 
7imp and the dotted line is a fit with the y/T function. 
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